Electron-Phonon Interactions in Polyacene Organic Transistors 
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We present a simple model for the electron-phonon interactions between the energy subbands in 
polyacene field-effect transistors and the vibrations of the crystal. We introduce a generalized Su- 
Schrieffer-Heeger model, arguing that the strongest electron-phonon interactions in these systems 
originate from the dependence of inter-molecule hopping amplitudes on collective molecular motion. 
We compute the electron-phonon spectral function a 2 F(lo) as a function of two-dimensional hole 
density and the coupling strength constant. Our results are in agreement with the sharp onset of 
superconductivity near half-filling discovered in recent experiments by Schon et al. and predict 
an increase of T c with pressure. We further speculate on the implications that the observation of 
the quantum Hall effect in these systems has on the effective band mass in the low carrier density 
regime. 
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Present studies of high-mobility organic transis- 
tors by Batlogg and collaborators have shown an im- 
mense richness of tunnability of their electronic proper- 
ties with gate voltage of in a single device. New possibil- 
ities for studying the physics of coherent band quasipar- 
ticle transport in organic semiconductors have been open 
by the high quality of these organic single crystals and 
of their interfaces with AIO3 dialectics. To date, how- 
ever, the analysis of organic field-effect-transistor elec- 
tronic systems has been hampered by the absence of sim- 
ple and reliable models for their electronic quasiparticles 
and for the interactions of these quasiparticles with each 
other and with vibrations of the host lattice. Indeed the 
discovery of coherent band quasiparticle properties and 
of the quantum Hall effect in these systems, which are 
usually thought of as being complex and relatively dis- 
ordered, has been one of the major surprises that has 
emerged from recent materials quality advances. 

At temperatures above « 2QK the transport electronic 
properties of polyacene semiconductors are well described 
by small polaron theory ||. However, this small polaron 
theory is based on localized molecular orbitals and is 
therefore unable to describe the low-temperature band 
quasiparticle behavior seen in Battlogg et al. 's samples. 
Schubnikov-de Haas transport measurements in these 
samples demonstrate a quasiparticle mean-free-paths in 
excess of 1000 lattice constants, demonstrating that the 
low-temperature regime can be described only by starting 
at the opposite limit and building a theory based on de- 
localized band quasiparticles states. Nevertheless, the in- 
teractions of quasiparticles with lattice vibrations remain 
strong, and are presumably responsible for superconduc- 
tivity. We propose a simple tight-binding model for the 
quasiparticle bands originating from the 7r-molecular or- 
bitals and for the interaction of these bands with the host 
crystal vibrational excitations. Our theory is parameter 
free once the band mass has been specified Q . 

We focus, for simplicity, on the polyacene compound 
anthracene, which is monoclinic (P2j/ ) with two ba- 
sis molecules whose orientations are related by a glid- 



ing plane symmetry. We start by considering the 
field-induced energy bands obtained neglecting electron- 
phonon interactions. The key question we need to ad- 
dress is the number of 2D subbands that are occupied 
at a particular 2D density, nio- We find that all carri- 
ers reside in a single subband up to much higher density 
in these transistors than in their inorganic counterparts, 
principally because of the relatively small dialectic con- 
stants (e ~ 3.5) and the large effective in plane band mass 
(toj_ = 1.5m e ) and in spite of the rather large c-direction 
effective band mass m z ~ 3m e Our conclusion is 
based partially on envelope-function density-functional 
calculations for which the local-density-approximation 
(LDA) can be used to estimate many-body effects Q 
that favor a single subband. In our LDA calculations 
we find that the second subband is first populated at 
n2D ~ 5 — 7 x 10 14 cm -2 , a density just larger than 
the highest achieved in experimental systems. Since the 
experimental signatures of second subband occupation 
would be unambiguous and no anomalies in 2D density 
dependence have been reported, we conclude that a sin- 
gle 2D subband is occupied up to the highest densities 
and that the lowest subband is strongly localized in the 
top layer for n,2D larger than i=s 10 13 cm -2 . 

Because Huckel model intra-molecular hopping pa- 
rameters have no linear dependence on low energy 
intra-molecular vibration normal coordinates, our model 
Hamiltonian includes only inter-molecular hopping pa- 
rameters H . Then to model the electron-phonon in- 
teraction, we expand each a — b plane hopping param- 
eter to first order in the twelve coordinates that describe 
rigid rotations and displacements of neighboring poly- 
acene molecules (as it is usually done in the Su-Schrieffer- 
Heeger model M): 
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where U^, m is the generalized displacement coordinate, 
fj, = 1,2 is the molecular basis index, m — 1,2,3 de- 
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note displacements of the molecular center of mass along 
x,y and z directions, and m = 4, 5, 6 denote angular 
displacements around the 1,2,3 principal axes angular 
displacements around the 1,2,3 principal axes of each 
molecule. The electron-phonon interaction parameters 
t^ m = dt/du^^m are calculated based on the assumption 
of proportionality between hopping integrals and overlap 
integrals between HOMO or LUMO orbitals on adjacent 
molecules. These are calculated using standard Hiickel 
approximation HOMO 7r orbitals and hence appropriate 
to the hole systems on which we focus. 

We next present expressions for the interactions terms 
in the Hamiltonian between phonons of the host crys- 
tal and the 2D band quasiparticles. The phonon fre- 
quencies and polarization vectors are obtained by solv- 
ing the standard secular equation involving the dynam- 
ical matrix. These eigenmodes are in general a mixture 
of displacements and rotations. We compute the dynam- 
ical matrix following the procedure outlined in Ref. 0. 
The phonon density-of-states and dispersion curves that 
emerge from these calculations have been reported else- 
where and are omitted here for the sake of brevity ||. 
The rigid-molecule vibration approximation used here is 
a convenient but inessential approximation that can be 
complemented if necessary by including coupling to im- 
portant isolated molecule vibrations; it is however reason- 
ably accurate for low-frequency vibrations of anthracene 
becoming less reliable for larger polyacenes. 

Combining these ingredients we finally write the 
Hamiltonian (ignoring the Coulomb interaction term) as 
^ = ^2D + W 2 Dc-vib + H V ib, where 
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and 
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In Eq. | f a .f}{q z ) = J2 Zl af*a^e l ^ z ' is the form factor of 
the subband (indexed by i, j) tight-binding wave function 
and the brackets indicate reduction to the 2D BZ. In 
our calculations we use the strict 2D limit (f a , a (Qz) — 
1) because of the strong quantum confinement found in 
LDA and Hartree calculations. 

We next define the electron-phonon interaction spec- 
tral function as 

E El^i;M(k,[k'-kU,*/)| 2 
x<5(e;(k) - /i)«(e i (k / ) - n)8(w - M[k' - k], fe ) (4) 



where g(fi) — 5Z k £(ej(k) — A 1 ) 1S the electronic den- 
sity of states, and [k' — k] denotes the projection of 
the Q = (k' — k, q z ) vector into the three dimen- 
sional BZ. The results of these calculations for densi- 



ties n 2 d = 4.0 x 10 14 cm~ 2 , n 2 u = 1-4 x 10 14 cm" 



and 



«2D = 6.7 x 10 13 cm -2 are shown in Fig. [l]. The peak 
locations in the low-frequency rigid molecule vibration 
regime agree with those observed in infrared absorption 
and tunneling data in pentacene crystals ||. The mass 
enhancement factor A, obtained from the integration of 
2Fa 2 (oj)/ui, close to half filling is 0.25. Its sharp decay 
away from half-filling (driven by the density of states) 
helps explain the sharp superconductivity onset observed 
in the experiments, and as a consequence we expect a 
suppresion of superconductivity beyond half-filling from 
this simple theory. We estimate the superconducting crit- 
ical temperature by the [dimension independent] BCS ex- 
pression fcsT c sw huz) exp[— 1/Al. Such expression, using 
huju/kB ~ 150K from Ref. M, yields ~ 3Kin agree- 
ment with the experiments Q. We note that the inclu- 
sion of the second-nearest-neighbour hopping tends to 
push the density of states peak to higher energies (den- 
sities). Adding a finite life time of the quasiparticles, 
will broaden these peaks, hence creating a plateau in the 
n2D-dependence of T c shifted to the right of half filling 
as observed in the experiments . This last effect would 
lower also the T c calculated at n — 4.00 x 10~ 12 cm~ 2 , 
worsening experimental agreement. However, given that 
Ep does not satisfy Ep >> Ku) max , such estimates of T c 
from A must be considered qualitative. 

In summary we have presented a theory of the low- 
temperature quasiparticle bands in polyacene field effect 
transistors, and of the electron-phonon coupling with the 
host molecular crystal. Our calculations indicate that 
the quasiparticles lie in a single 2D tight-binding band 
up to the highest densities that have been achieved ex- 
perimentally at present and that the most important 
electron-phonon interactions arise from the influence of 
approximately rigid molecular translations and rotations 
on hopping between HOMO and LUMO orbitals on ad- 
jacent molecules. This picture implies a dependence of 
superconducting critical temperature on molecular lat- 
tice constant which contrasts with the case of doped 
fullerene superconductors. For the fullerenes, the impor- 
tant electron-phonon interactions are intra-molecular so 
that T c depends on inter-molecular hopping only through 
the density-of-states. Decreasing the lattice constant, 
increases hopping, and hence decreases the density-of- 
states and T c ||. Here decreasing the lattice constant 
will also strengthen the important electron-phonon inter- 
actions. According to our theory, the latter effect dom- 
inates and T c will increase with decreasing lattice con- 
stant. We find that phonon-mediated electron-electron 
interactions in polyacene molecular crystals are strong 
only when the 2D tight-binding band is close to half- 
filling and its density of states is relatively large. This oc- 
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curs for 2D densities comparable to 4 x 10 14 cm -2 , where 
superconductivity turns on relatively abruptly. At lower 
carrier densities, A <C 1 and electron-electron interac- 
tions are dominated by repulsive Coulomb interactions, 
consistent with the occurrence of the fractional quantum 
Hall effect. We also note that the simple fact that the 
fractional quantum Hall effect is observed points to the 
distinct possibility that the band mass can be highly re- 
duced at those lower carrier densities since, otherwise, 
from rough estimates of the known experimental values 
at higher carrier densities, the r s value would be too large 
to allow anything but a Wigner crystal ground state. 
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FIG. 1. a 2 F{io) for n 2d = 4.0 x 10 14 cm" 2 (solid line), 
n 2d = 1.4xl0 14 cm -2 (dotted line), and n 2d =6.7xl0 13 cm -2 
(dashed line). 
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